ASYMPTOTIC VALUATIONS OF SEQUENCES 
SATISFYING FIRST ORDER RECURRENCES 
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Abstract. Let t n be a sequence that satisfies a first order homogeneous 
recurrence t n = Q(n)t n -i, where Q £ Z[n]. The asymptotic behavior 
of the p-adic valuation of t„ is described under the assumption that all 
the roots of Q in Z/pZ have nonvanishing derivative. 



1. Introduction 
The p-adic valuation v p (x), for x € Q, x ^ 0, is denned by 
(1.1) x = p v 



b 1 

where a, b G Z and p divides neither a nor b. The value f p (0) is left unde- 
fined. 

In this paper we establish the asymptotic behavior of the p-adic valuation 
of sequences that satisfy first order recurrences 

(1.2) t n = Q(n)t n . u n>l, 

where Q is a polynomial with integer coefficients. Among all the positive 
integer zeros of Q, let v be the maximum modulus. Take uq > v. Then the 
recurrence (|1 .2j) is started at this index no- This ensures the non- vanishing 
of t n . Without loss of generality, we always assume no = and to = 1. 
We also adopt the notation t n (Q) while refering to the sequence defined by 

<d2D. 

The identity 

n 

(1.3) "pWO)) = E*))' 

i=l 

shows that only the zeros of Q in Z/pZ contribute to the value of v p {t n {Q)). 
The main tool of our asymptotic analysis will be Hensel's lemma. The 
version stated here is reproduced from [3]: 

Lemma 1.1 (Hensel). Let f{x) G Z p [x] be a polynomial with coefficients 
in the p-adic integers Z p . Write f'{x) for its formal derivative. If f{x) = 
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mod p has a solution at satisfying f'{a\) ^ mod p, then there is a unique 
p-adic integer a such that f(a) = and a = a\ mod p. 

We now state our main theorem. This result is an asymptotic description 
of the valuation of the sequence t n , defined by (jl.2[) . 



Theorem 1.2. Let Q € Z[n]. Assume each of the roots of Q satisfies the 
hypothesis of Hensel's lemma. Let z p denote the number of roots of Q in 
TLjpTL, that is, 

(1.4) z p := \{b e {1, 2, . . . ,p} : Q(b) = mod p}\. 
Then the sequence {t n }, defined in 111. Sty , obeys the estimate 

Z Ti 

(1.5) v p {t n ) ~ — — j- as n — » oo. 

Motivation. The most elementary example is Q(x) = x. Theorem 11.21 
yields v p {n\) ~ n/(p — 1). This follows from the classical formula of Legendre 

n s f n n - sjn) 

(1.6) f p (n!) = 



p—1 

where s p (n) is the sum of the digits of n in base p. 

Our motivation for Theorem 11.21 comes from the study of the sequence 
{x n } defined by 

n 

(1.7) x n = tan tan -1 k, n > 1. 

fc=i 

This same sequence satisfies the recursive relation 

x n _i + n 



(1.8) x r 



1 - nx n -\ 



with initial condition x\ = 1. The first few values are {1, —3, 0, 4, — 
and in [2] it was conjectured that x n ^ for n > 4. Later this was proved 
in [1] using the 2-adic valuation of x n . The sequence {x n } was linked in pQ 
to 

(1.9) w n := (1 + 1 2 )(1 + 2 2 )(1 + 3 2 ) • • • (1 + n 2 ), 
which can be condensated as 

(1.10) w n = (l + n 2 K-i- 

This corresponds to Q(x) = x 2 + 1 and it fits into the type of recurrences 
considered here. 

Section [2] contains the proof of Theorem 11.21 and Section [3] presentes ex- 
amples illustrating the main result. In the last section we propose some 
future directions. 
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2. The proof 

In the proof we assume that Q has no roots in N U {0}. The general 
situation can be reduced to this one by a shift of the independent variable. 

The conclusion of Theorem 11.21 is trivial if z p = 0, so we assume z p > 0. 
Denote by b%, b%,--- ,b Zp the zeros of Q in Z/pZ. The definition of t n yields 

n 

(2.1) v p {t n ) =5>p(Q(*)). 

8=1 

All sums below are assumed to run from i = 1 to n. 

Only the indices congruent to bj modulo p contribute to (|2.1j) . thus 

(2.2) u p (t n )= Yl "?(<?(»)) + -+ E 

i=6l mod p i=t>z p mod p 

where 1 < i < n. For fixed j G {1, 2, . . . , z p }, we consider the term 

(2.3) Yl MQ(i))- 

i=bj mod p 

Hensel's lemma produces a p-adic integer 

(2.4) = p jt0 + (3 jA p + ■■■ + (3 jjkP k + ■■■ 

such that (3j t k € {0, 1, • • • , p — 1}, /J^q = 6j mod p and Q((3j) = 0. Observe 
that if the representation (|2.4p were finite, then /3j would be a non-negative 
integer root of Q. This possibility has been excluded. Introduce the notation 

(2-5) 7 i>a := + p/3j-,i + p 2 /? ii2 + ■ ■ • + p s (3j, s p s - 



Definition 2.1. For n € N, Zet 

(2.6) r n = Max {j : p 3 divides some Q{i) for 1 < i < n}. 

Lemma 2.2. TTie sequence r n — > oo as n — » oo. Moreover, for large n, we 
have p Tn < n deg ^ +l , hence r n = O(logn). 

Proof. Hensel's lemma shows that 7j jS satisfies Q(7j, s ) = Omod p s+1 . For 
any given M > 0, choose an integer s > M. Taking n > 7j,s— i we have that 
* := 7?>-i ^ {!) 2, ••• ,n} and The definition of r n implies that 

r n > s > M. Therefore r n — > oo as n — > oo. Now observe that p Tn divides 
|Q(i)| for some 1 < i < n. The estimate 

(2.7) p rn < \Q(i)\ < Max{|Q(l)|, • • • , \Q(n)\} < Cn deg ^ 

gives the upper bound on r n . The constant C depends only on the coeffi- 
cients of Q. □ 

Now 

E v p(q(}))= E 1+ E ! + •••+ E x > 

!E6j mod p i=7j i omodp »=7j,l mod p 2 *=7j,r n -i mod p r ™ 
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where all sums range over 1 < i < n. The bound 



(2.8) 
yields 



n 
p s 



^ £ is 

i=fj,s mod p 



n 



p s 



+ 1 



E ^(3(0) > 



i=bj mod p 



1 1 



7^ 



1 



pl'n 



+ — + •■■+ „ 



n 



p — 1 



1-P 



Therefore 



p-1 



(p - iK, 



i=bj mod p 



and passing to the limit we conclude that 



(2.9) 



lim inf 



p — 1 



n— »oo n 



E > i- 



i=bj mod p 



Similarly, using the upper bound in (|2.8p we obtain 

i=bj mod p 

and it follows that 

p — 1 



P — 1 



(2.10) 



lim sup ; 



n— >oo Tl 



E ^(o(0) < i- 



i=bj mod p 



Therefore, Theorem 1 1 . 2 1 has been established. 

3. Examples 

In this section we present some examples illustrating Theorem 11.21 

Definition 3.1. Given a polynomial Q and a prime p, we say that a G TLjpTL 
is a Hensel zero of Q if Q(a) = mod p and Q'(a) ^ mod p. The prime p 
is called a Hensel prime for Q if all the zeros of Q in Z/pZ are Hensel zeros. 
We also require that Q has at least one zero in Z/pZ. The asymptotic zero 
number is defined (provided it exists) by the limit 

(3.1) N P (Q) := lim {p - . 

Theorem 11.21 is restated as follows: 



Theorem 3.2. Let p be a Hensel prime for Q. Then N p (Q) = z p . 
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Note. The examples will show pairs (Q,p) for which N p {Q) ^ N. An 
appropriate interpretation of this number is lacking in these cases. 
In the examples described below we present the normalized error 

(3.2) err p (n; Q) := z p n - (p - l)v p (t n {Q)) 
and the relative error. 

(3.3) relerr p (n; Q) := err p (n; Q) — err p (n — 1; Q). 

Certain regular structure of this function, as seen in Figure [3j will be ana- 
lyzed in a future report. 

Example 1. Let Q(x) = x 5 + 2x 3 + 3. Then p = 5 is a Hensel prime 
for Q. Indeed, the only zeros of Q in Z/5Z are a = 3 and a = 4 and 
Q'(a) ^ mod 5. Theorem 11.21 gives 

n 

(3-4) vsMQ)) ~ 2" 

Figure [1] shows the valuation u p (t n (Q)). Figure [2] and [3] depict patterns in 
the normal and relative error, respectively. 




Figure 1. The valuation u 5 (t n ) for Q(x) = x 5 + 2x 3 + 3. 



Example 2. A direct calculation shows that, among the first 20000 primes, 
p = 3, 11 and 29 are the only non-Hensel primes for Q(x) = x 5 + 2x 3 + 3. We 
now describe the asymptotic behavior of f p (t n (Q)) in each of these cases. 
The polynomial Q factors as 



(3.5) 
where 

(3.6) H(x) 
and the valuation splits as 



x 5 + 2x 3 + 3 



(sc + l)H(x) 



x 3 + 3x 2 - 3x + 3 



(3.7) 



v p {t n {Q)) = v p (t n (x + 1)) + u p {t n (H{x)). 
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-5 - 




Figure 2. The normalized error when p = 5 and Q(x) = 
x 5 + 2x 3 + 3. 




Figure 3. The relative error when p = 5 and Q(x) = x 5 + 
2x 3 + 3. 

Theorem 11.21 gives v p (t n (x + 1)) ~ n/{p — 1), so it remains to evaluate 
u p (t n (H)). 

The prime p = 3. In this case and 1 are zeros of H in Z/3Z, and only 1 
is a Hensel zero. Observe that 

(3.8) v 3 (t n (H))= E ^(#(j))+ £ Mm)- 

j=0 mod 3 j=l mod 3 

Since 1 is a Hensel zero, the argument in the proof of Theorem 11.21 implies 
that 

(3.9) Yl Mm)~^- 

j=l mod 3 

To analyze the first sum in (|3.8p . note that 

(3.10) H(3k) = 81/c 4 - 27k 3 + 27k 2 -9k + 3. 
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Thus, u 3 (H(3k)) = 1 for k G N. We obtain that 

n 
3' 



(3-11) £ r/ 3 (i?(i)) " 



j=0 mod 3 

and then u 3 (t n (H)) ~ Therefore i*(t„(Q)) ~ ^ and JV 3 (Q) = |. 

The prime p = 11. For this prime, although Theorem 11.21 does not apply 
to Q itself, it is applicable to both factors x + 1 and H{x). And, we deduce 

,^nn n In 3n 

(3.12) vn(tJQ)) H = — • 

Therefore N n (Q) = 3. 

The prime p = 29. In order to find the asymptotic behavior of V2§{t n {H)), 
observe that 14 is the only zero of H in Z/29Z and 

(3.13) #(29A; + 14) = 36221 + 303601 k + 956217A; 2 + 1341395fc 3 + 707281/c 4 . 

The valuations of the coefficients in H(29k + 14) are 1, 2, 2, 3, and 4, re- 
spectively. Therefore V2${H{j)) = 1 if j = 1 mod 29 and otherwise. We 
conclude that 

Tl 

(3-14) M*n{H)) ~ ^7^- 

Therefore i/ 29 (<n(Q)) ~ If and iV 29 (Q) = §. 
Example 3. The polynomial 

(3.15) Q(x) = x 8 + x 5 + x 3 + 1 = (x 3 + l)(x 5 + 1) 
does not have a Hensel prime. This follows from 

(3.16) gcd(Q0r),Q'(x)) = x + 1, 

so that, for any prime p, we have that p — 1 is a zero of Q in Z/pZ and 
Q'{p — 1) = 0. Naturally we have 

(3.17) Vp{tn{Q)) = Vp(tn(x 3 + 1)) + lS p {t n (x 5 + 1)). 

The asymptotic behavior of v p (t n (Q)) is discussed next. 
Lemma 3.3. Let p be an odd prime and x ^ 1. Then 

{0 if x ^ 1 mod p 

1 + f p (x — 1) if x = 1 mod p. 



(3.18) z/ p (x p - 1) 



Proof. The first part is clear from the congruence x p = x mod p. To verify 
the second assertion, write x = kp + 1 and observe that 

p 

(3.19) u p {x p - 1) = i/ p ' " 

For r > 1, the p-adic valuation of each term in the sum is greater than 
2 + v p (k). When r = 1, it is exactly 2 + v p (k). Then, putting k = ^y- 
verifies the assertion. □ 
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Corollary 3.4. Let p be an odd prime and x 6 Z, a; 7^ 1. Define 

(3.20) T p (x) = x p ~ l + x p ~ 2 + ■■■ + !. 

Then 



(3.21) v P (T p (x)) 



if x 1 mod p 

1 if x = 1 mod p. 



Corollary 3.5. Lei p be a prime and idZ, TTien 

(3.22) ^+d={" . _ 

I 1 + v p (x + 1) x = — 1 mod p. 

Proof. Replace x by — x in Lemma 13.31 □ 
Corollary 3.6. Let p be a prime and x £ Z, x 7^ — 1. Define 

(3.23) 5 p (x) = x p - x - x p " 2 + x + 1. 

T/ien 



(3.24) fp(S p (x)) 



if x ^ —1 mod p 

1 if x = —1 mod p. 



The number of roots of x 9 + 1 = mod p, that is, z p (x 9 + 1) stated in the 
Lemma below appears at the end of Section 8.1 of [3]. 

Lemma 3.7. Let p and q be primes. The number of solutions of the con- 
gruence x p + 1 = mod q is gcd(p, q — 1). 

Corollary 3.8. Let p be an odd prime. Then 

(3.25) i/ p (t w (gP±l))~ ^'^ly 
If q is a prime, 9 ^ p, then 

(3.26) uSn^±l))- 9Cd{P ^\ 1)n . 

Proof. Theorem 11.21 gives v p (t n (x + l)) ~ The expression for v p (S p (x)) 
yields z/ p (S p (x)) ~ n/p. The asymptotic behavior of v g (t n (x p ± 1)) follow 
directly from Theorem 11.21 □ 

We now complete the analysis of 

(3.27) v p {t n {Q)) = v p (t n {x 3 + 1)) + v p {t n (x 5 + 1)). 
If p 7^ 3 is a prime, then 



(3.28) ^ p (t n (x 3 + 1)) ~ Zp(x3+ 1 1)n 

p — 1 
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Similarly, for p 7^ 5 prime, we have 



(3.29) u p (t n (x 5 + 1)) ~ Zp(x + 1)n . 

p — 1 

Thus, (|3T25l) and (|339l) yield 



~ ^3(tn(x 3 + 1)) + V 3 (t n (x 5 + 1)) = ^ + ~ = ^. 

2 3 

Similarly, u^{t n {Q)) ~ 7n/10. 

Now let p 7^ 3, 5 be a prime. Theorem 11.21 now applies directly to give 

(3.30) MUQ)) ~ Mf!±ll±^f!±^ . 

p — 1 

Lemma 13.71 yields 

/oqi\ /+ /^\\ [gcd(3,p - 1) + gcd(5,p - 1)] n 

(3.31) Vp\pn{Q)) ~ i • 

p — 1 

The asymptotic zero number is given by 



(3.32) N p ((x 3 + l)(x 5 + l)) = < 



§ if p = 3 

f ifp = 5 

gcd(3,p- 1) +gcd(5,p- 1) if p ^3, 5. 



Example 4. Let p be an arbitrary prime and define 
(3.33) Ap[x) = (px + l) 2 ((p+l)x + l). 

A direct calculation shows that p is the only Hensel prime for A p . Therefore 

n 



(3.34) v p {t n {A p )) 

p — 1 

To compute the asymptotics for a prime q 7^ p, let Qi(cc) = px + 1 and 
Q2(^) = (p + l)ic + 1, and observe that 

(3.35) v q {t n {A p )) = 2v q (t n (Qx)) + v q {t n {Q 2 )). 

Theorem 1 1 . 2 1 applies to both Q\ and Q 2 - The case for Q\ is immediate since 
px + 1 = mod q has a unique solution. To evaluate v q {t n {Q2)) observe 
that the number of solutions of (p + l)x + 1 = mod q is or 1, according 
to whether q divides p + 1 or not. Thus 

3.36 Mtn(Ap)) ~ 1 p ' g; 

9 — 1 

where 

J 1 if g divides p + 1, 
I otherwise. 
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We conclude that 
(3.37) N q (A P ) = 




if p = q, 
P,q ifp^q- 



4. Future directions 

In this section we outline certain generalizations of the main result of the 
paper. 

A natural extension of Theorem 11.21 deals with the situation in which 
there is an element b € Z/pZ such that 

(4.1) Q(b) = Q'{b) = ■■■ = Q^- 1 ) (b) = mod p. 

The question of how the multiplicities of the roots enter in the asymptotic 
behavior of v p {t n {Q)) appears to be a salient quest, and this will be addressed 
elsewhere. 

Another interesting continuation of the ideas presented in this paper 
would be the study of p-adic valuation of sequences satisfying second or- 
der recurrences 

(4.2) t n = Qi(n)t„_i + Q 2 (n)*n-2, 

with polynomials Q\ and Q2- This problem includes, classically, the case of 
Fibonacci and Stirling numbers. 

Acknowledgments. The work of the third author was partially funded by 
NSF-DMS 0409968. The second author was partially supported as a grad- 
uate student by the same grant. 



References 

[1] T. Amdeberhan, L. Medina, and V. Moll. Arithmetical properties of a sequence arising 
from an arctangent sum. To appear in Journal of Number Theory, 2007. 

[2] G. Boros and V. Moll. Sums of arctangents and some formulas of Ramanujan. Scientia, 
11:13-24, 2005. 

[3] K. Ireland and M. Rosen. A classical introduction to Number Theory. Springer Verlag, 
2nd edition, 1990. 

[4] M. Ram Murty. Introduction to p-adic Analytic Number Theory, volume 27 of Studies 
in Advanced Mathematics. American Mathematical Society, 1st edition, 2002. 

Department of Mathematics, Tulane University, New Orleans, LA 70118 
E-mail address: tamdeber@tulane.edu 

Department of Mathematics, Tulane University, New Orleans, LA 70118 
E-mail address: lmedina@math.tulane.edu 



Department of Mathematics, Tulane University, New Orleans, LA 70118 
E-mail address: vhm@math.tulane.edu 



